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Frequency-domain  methods  are  used  to  extract  the  open-loop  dynamics  of  the  XV-IS  tilt- 
rotor  aircraft  from  flight  teat  data  for  the  cruise  condition  (V  -  170  knots) .  The  frequency 
responses  are  numerically  fitted  with  transfer-function  forms  to  identify  equivalent  modal 
characteristics.  The  associated  handling  quality  parameters  meet  or  exceed  Level  II,  Cate¬ 
gory  A,  requirements  for  fixed-wing  military  aircraft.  Step  response  matching  is  used  to 
verify  the  time-domain  fidelity  of  the  transfer-function  models  for  the  cruise  and  hover  flight 
condltioiui.  The  transient  responses  of  the  model  and  aircraft  are  in  close  agreement  in  all 
cases,  except  for  the  normal  acceleration  response  to  elevator  deflection  in  cruise.  This  dis¬ 
crepancy  is  probably  due  to  the  unmodeled  rotor  rpm  dynamics.  The  utility  of  the  frequency- 
domain  approach  for  dynamics  Identification  and  analysis  is  clearly  demonstrated. 

1.  Introduction 

The  identification  of  XV-15  tilt-rotor  dynamics  (Fig.  1)  from  flight  test  data  is  an 
extensive  ongoing  effort  to  support  the  development  of  the  next  generation  of  tilt  rotors— the 
Joint  services  V/STOL  aircraft  (JVX) . 

The  key  concerns  of  the  effort  are  the 
documentation  of  open-loop  XV-15 
dynamics,  and  the  validation  of 
generic  tilt-rotor  models  (Refs.  [1] 
and  [2]).  A  frequency  domain-based 
Identification  approach  was  developed 
and  successfully  applied  for  the  hover 
flight  condition  (Ref.  [3]).  Transfer 
function  models  describing  the  open- 
loop  response  characteristics  of  the 
XV-IS  aircraft  were  extracted  and  com¬ 
pared  with  the  simulation  characteris¬ 
tics.  Reference  [3]  presents  a 
detailed  description  of  the  frequency- 
domain  methodology  and  the  results  for 
the  hover  flight  condition. 

As  in  all  identification 
research,  a  key  concern  of  this  effort 
is  the  fidelity  of  the  extracted 
models  for  input  forms  other  than 
those  used  in  the  identification  pro¬ 
cess.  The  flight  testing  technique  of 
Ref.  [3]  uses  a  pilot-generated  swept 
sine  wave  input  to  excite  the  vehicle 
dynamics.  This  yields  an  excellent 
Identification  of  frequency  responses 
and  transfer  functions  for  sinusoidal- 
like  inputs.  Such  information  is 
useful  for  frequency  domain-based 
handling  quality  specifications  such 
as  the  Mil  Handbook  for  military  air¬ 
craft  (Ref.  [4]).  However,  for  time- 
domain  specifications  such  as 
Mil-H-8501  (Ref.  [5])  and  Mll-F-83300 
(Ref.  [6]),  criteria  are  largely  based 
on  responses  to  step  inputs,  so 
extracted  models  must  accurately 
reflect  these  characteristics  as  well. 

Figure  1.  The  XV-15  Tilt-Rotor  Aircraft,  (a)  Cruise 
configuration,  (b)  Hover  configuration. 


Thl*  paper  revlawa  the  frequency-douln  based  nethodology  and  discusses  the  identifica¬ 
tion  of  transfer  function  aodels  for  the  cruise  flight  condition  (170  knots).  This  is  a  good 
llBitlng  case  for  co^>arlson  with  the  previous  hover  results  (Ref.  [3]).  TlM-douln  astchlns 
is  presented  to  verify  the  step  response  characteristics  of  the  extracted  aodels  for  the  hover 
and  cruise  fll^t  conditions.  These  results  show  the  utility  of  relatively  slnple  transfer 
function  aodels  for  handling  qualities  and  control  system  applications. 


2.  Review  of  Identification  Methodology 

This  section  reviews  the  frequency-domain  Identification  and  time-domain  verification 
techniques.  The  details  of  the  Identification  approach  are  extensively  discussed  in  Ref.  (31 
and  are  only  outlined  in  this  paper.  *  ‘  ^ 

Frequency-domain  identification  is  based'on  the  spectral  analysis  of  input  and  output 
time  histories  using  Fast  Fourier  Transform  techniques.  This  analysis  produces  describing 
functions  which  are  (coiiq>lex-valued)  linear  descriptions  of  the  Input-to-output  processes.  The 
identification  results  are  presented  in  Bode  plot  format:  magnitude  and  phase  versus  frequency. 
System  bandwidth  and  effective  time  delay.  Important  metrics  in  current  handling  qualities 
specifications,  can  be  read  directly  from  these  plots.  Tabulated  frequency  response  results 
are  fitted  with  analytical  transfer  function  forms  to  extract  modal  characteristics  for  handling 
qualities  specifications  given  in  terms  of  lower  order  (equivalent)  system  models.  Also,  the 
trsnsfer-functlon  models  can  be  driven  with  step  inputs  to  extract  familiar  time  domain  metrics 
such  as  rise  tiam,  overshoot,  and  settling  time. 


The  swept  sine  wave  (frequency-sweep)  input  is  a  good  excitation  for  the  frequency-domain 
identification  approach.  This  excitation  results  in  bounded  and  reasonable  excursions  of  the 
aircraft,  suitably  exciting  the  Important  rigid  body  modes  over  the  entire  frequency  range  of 
interest  (0.2-d.0  rad/sec).  The  input  is  generated  by  the  pilot  in  one  axis  at  a  time,  with 
minimal  regulation  of  the  remaining  axes.  Starting  from  a  trim  condition,  the  pilot  first  exe¬ 
cutes  two  20  sec  period  inputs  to  ensure  good  low-frequency  identification.  Then  the  frequency 
of  the  Inputs  is  slowly  increased  up  to  a  maximum  of  about  6  rad/sec,  yielding  a  total  run 
length  of  about  90  sec.  Repeat  runs  are  executed  to  allow  concatenation  of  the  90  sec  records 
in  the  analysis;  this  technique  minimizes  the  effects  of  random  noise. 


When  the  bare  airframe  characteristics  of  the  vehicle  are  highly  unstable,  as  in  the 
XV-15  in  hover,  execution  of  a  90  sec  frequency  sweep  on  the  open-loop  vehicle  is  clearly  not 
practical.  Subject  to  some  important  qualifications,  the  bare  airframe  transfer  functions  can 
be  extracted  from  the  closed-loop  flight  data.  Referring  to  Tig.  2,  identification  of  the  open- 

loop  transfer  function  (q/d,)  requires 

..  (BARE-  the  measurement  of  the  output  (q)  and 

(PILOT)  AIRFRAME)  (OUTPUT)  surface  deflection  (6^)  signals.  This 

(INPUT)  i  - 1  total  surface  deflection  is  made  up  of 

q  components  from  the  pilot  (d^oj))  and  the 
stability  and  control  augmentation  sys¬ 
tem  (SCAS)  feedback  (df^) .  A  key 

requirement  is  that  the  total  surface 
deflection  must  contain  a  significant 
component  from  the  pilot  which  is  uncor¬ 
related  with  the  output.  The  use  of  the 
"prograoMd"  swept  sine  wave  input 
avoids  the  otherwise  troublesoise  result 
of  identifying  the  inverse  feedback 


(SCAS) 


transfer  function  (-1/Cp^)  from  passive 
tracking  tasks. 


When  the  vehicle  dynamics  are 

Figure  2.  Single  degree-of-freedom  pitch  response  highly  coupled,  multiple  control  deflec- 
*“***^*  tlons  will  occur  in  all  degrees  of  free- 

.  ^  ^  If  the  surface  deflections  are  not 

jUlly  oovrolated  with  each  other,  the  desired  single-input  transfer  function  can  be  isolated 
using  the  multi-input  multi-output  methodology  described  in  Ref.  [3).  This  restriction  was 
satisfied  for  the  hover  flight  condition  since  the  cross-control  inputs  supplied  by  the  pilot 
l^tgely  uncorrelated  with  the  primary  frequency  sweep  inputs.  Therefore,  it  was  possible 
to  Isolate  the  yaw  response  to  ailerons—  an  Important  source  of  coupling  for  the  hover  flight 
condition. 

f***  cruise  flight  condition,  open-loop  longitudinal  transfer  functions  were  also 
Idmtlflcd  from  the  closed-loop  tests  since  the  elevator  Inputs  supplied  by  the  pilot  were  uncor- 
relatsd  with  the  outputs  and  the  other  surface  deflections.  In  the  lateral  axis,  the  large 
"*®*’**  Inherent  closed-loop  stability  allowed  the  pilot  to  execute  frequency  sweeps  of  the 
slleron  without  applying  rudder  corrections.  However,  the  stability  augmentation  system  supplies 


rudder  inputs  which  are  highly  correlated  to  the  aileron  Inputs  and  roll  rate  outputs.  There¬ 
fore.  extraction  of  the  single-input,  single-output  open-loop  transfer  functions  from  closed- 
loop  flight  data  %ras  not  possible  for  this  flight  condition;  the  lateral-directional  tests  were 
repeated  with  the  SCAS  disengaged.  When  the  bare  airframe  dynamics  are  too  highly  coupled  or 
unstable  to  make  this  practical,  the  open-loop  dynamics  can  usually  be  extracted  from  the 
closed-loop  test  data,  as  was  successfully  accomplished  for  the  XV-IS  in  hover  (Ref.  [3]).. 

Once  the  flight  data  have  been  digitized  and  preprocessed.  Che  input,  output,  and  cross 
spectra  (Cxx(^) •  6yy(f) .  and  Gxy(f ) .  respectively)  are  calculated  using  modern  Chirp 
s-transform  methods  (Ref.  [3]).  Specific  transfer  functions  G(f)  are  obtained  from  the  ratio 
of  the  appropriate  cross  and  input  auto  spectra: 

-  ■  -  G  (f) 

G(f>  -  (1) 

XX 

and  are  presented  in  Bode  format.  The  coherence  function  (yt„)  defined  as: 


'xy  G  G 

'  XX ' '  yy ' 

is  a  good  indication  of  the  Input-Co-outpuC  linearity.  This  frequency-dependent  parameter  may 
be  Interpreted  as  that  fraction  of  the  output  spectrum  «diich  can  be  accounted  for  by  linear 
relation  with  the  input  spectrum.  When  the  process  under  investigation  is  perfectly  linear  and 
Che  spectral  estimates  are  noise  free,  the  coherence  function  will  be  unity  for  all  frequencies 
in  the  excited  input  spectrum  range.  A  value  of  the  coherence  function  less  than  unity  will 
result  from  nonlinearities  in  the  system,  input/output  noise,  or  cross-coupled  control  Inputs. 
The  magnitude  and  phase  responses  are  then  fitted  with  analytical  transfer  function  models  Co 
obtain  closed  form  descriptions  of  Che  input-to-output  processes.  In  order  to  obtain  a  unique 
fit  of  the  frequency  responses,  certain  physical  restraints  on  comnonality  of  the  transfer 
function  denominator  factors  are  Imposed. 

Current  military  handling  quality  specifications  (Ref.  [4])  limit  the  allowable  mismatch 
between  the  transfer  function  fits  and  the  frequency  response  data.  However,  even  small  decibel 
deviations  in  the  frequency-domain  can  produce  surprisingly  large  discrepancies  in  time-domain 
correlations.  Time-domain  verification  first  requires  converting  the  transfer  functions  to 
canonical  form: 


x(t)  =  Fk(t) 
2(t)  -  Cx(t) 


+  Gu(t)^ 


-a,  -aj  -a^ 


C  -  [»>,  bj  •••  b„  0  ...  O] 


and  a,,  aj^.  .  .  ..  a^_j  are  the  denominator  (ascending)  coefficients  and  b^.  b^ . 

are  the  numerator  (ascending)  coefficients  of  the  transfer  function  model  being  evaluated. 
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This  model  Is  driven  with  step- Input  data  from  the  flight  tapes,  and  the  resulting 
responses  are  then  compared  to  the  aircraft  data  for  tliae-domaln  fidelity.  This  approach 
allows  the  verification  of  transfer  functions  associated  with  the  atdble  degrees  of  freedom, 
since  spurious  model  inputs  (due  to  inexact  initial  condition  matching  and  turbulence  regula¬ 
tion)  do  not  cause  a  divergence  of  the  transients  with  time.  However,  for  the  wnatable  degrees 
of  freedom,  these  spurious  inputs  cause  a  rapid  divergence  of  the  flight  and  model  responses. 
For  such  cases,  conclusive  quantitative  statements  on  the  step  response  fidelity  of  the  trans¬ 
fer  function  models  cannot  be  made,  even  though  the  qualitative  form  of  the  response  is  correct 
and  small  differences  in  the  predicted  level  of  instability  are  probably  not  Important. 


3.  Dynamics  Identification  for  the  Cruise  Flight  Condition 

In  this  section,  the  characteristics  for  the  cruise  flight  conditio.^  are  considered: 

V  -  170  knots  (indicated),  nacelle  incidence  >  0*,  altitude  •  8000  ft.  The  longitudinal  and 
lateral  flight  dynamics  are  fully  decoupled  in  this  condition.  The  primary  longitudinal  bare- 
alrframe  transfer  functions  of  interest  are  pitch  rate  and  norsml  acceleration  responses  to 
elevator  deflection,  q/j^  and  respectively.  The  liportant  lateral-directional  transfer 

functions  are  roll  rate  response  to  aileron,  p/S^,  and  sideslip  (at  the  c.g.)  response  to 
rudder,  Bgg/«r. 

The  following  discussion  of  longitudinal  dynamics  identification  and  analysis  is  pre¬ 
sented  in  detail  to  illustrate  the  frequency-domain  methodology  and  to  compare  the  longitudinal 
results  for  the  hover  flight  condition  presented  in  Ref.  [3]. 


Longitudinal  Dynamics 

Since  the  longitudinal  and  lateral  dynamics  are  fully  decoupled  and  elevator  surface 
deflections  are  measured  directly,  the  flight  testing  for  the  longitudinal  axis  was  conducted 
with  the  SCAS  engaged.  This  was  not  considered  essential,  homver,  since  the  dynamics  for  this 
flight  condition  are  very  stable. 


The  longitudinal  stick  displacements  for  two  concatenated  frequency  sweeps  are  shorn  in 
Fig.  3.  As  previously  discussed,  each  run  lasts  roughly  90  sec  and  is  initiated  with  two  20-sec 

cycles  to  ensure  good  low-frequency  iden- 
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3.  Two  longitudinal  stick  frequency  sweeps 
in  cruise. 


Figure  4.  Elevator  surface  input  autospectrum 
(G6e6,) • 


tiflcation.  The  sinusoidal  stick  deflec¬ 
tion  is  very  regular  with  a  nearly  con¬ 
stant  amplitude  of  roughly  ±5Z.  The 
elevator  signal  is  also  very  regular  with 
a  nearly  constant  aiq>litude  of  2  deg. 

These  data  are  very  typical  of  the 
transition-flight  and  forward-flight 
condition  results.  The  relative  stabil¬ 
ity  of  the  longitudinal  and  lateral  axes 
made  execution  of  the  frequency  sweeps  a 
very  rapid  and  acceptable  flight-test 
procedure.  The  input  autospectrum  for 
the  concatenated  elevator  time  histories 
is  shown  in  Fig.  4.  This  spectrum  is 
fairly  flat  at  low  frequencies,  dropping 
off  slightly  for  the  higher  frequency 
inputs.  For  frequencies  outside  the 
plotted  range,  the  magnitude  drop-off  is 
much  more  pronounced. 

The  pitch-rate  signal  (conditioned 
with  a  2.5  Hz  low-pass  filter  for  output 
presentation  only)  is  shown  in  Fig.  5. 

The  pitch-rate  excursions  are  very  regu¬ 
lar  with  a  roughly  constant  peak-to-peak 
amplitude  of  5  deg/sec.  These  oscilla¬ 
tions  are  considered  very  acceptable  to 
the  research  pilots,  and  are  comparable 
with  those  encountered  during  the  hover 
flight  test.  The  associated  output  auto¬ 
spectrum  for  the  concatenated  pitch-rate 
time  histories  are  shown  in  Fig.  6.  The 
output  spectrum  follows  the  input  spectrum 
at  low  input  frequencies  and  exhibit  the 
familiar  20  dB/decade  roll-off  at  higher 
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frequencies  (which  Is  typical  of  rigid' 
body  motion) .  The  cross-spectrum 
between  elevator  surface  and  pitch 
rate  Is  shown  In  Fig .  7  and  demon¬ 
strates  a  similar  characteristic. 

The  elevator  surface-to-pitch 
rate  transfer  function,  q/d^.  Is  shown 
in  Figs.  8(a)  and  8(b).  The  smooth 
spectral  data  over  the  majority  of  the 
frequency  range  are  typical  of  the 
excellent  results  obtained  from  the 
Chirp  ■z-t-ranefcrm- algorithm.-  Due- to  ■ 
the  short  period-mode  excitation, 
there  Is  a  peak  In  the  pitch-rate 
response  at  about  2.0  rad/sec,  with 
an  associated  phase  lag  of  45  deg.  At 
higher  frequencies,  the  magnitude 
response  rolls  off  at  20  dB/decade  and 
the  phase  shift  approaches  -90  deg 
owing  to  the  K/s  pitch-rate  charac¬ 
teristic  and  negligible  flexibility/ 
servo-lag  effects.  The  drop  In  the 
magnitude  response  and  associated 
positive  phase  response  for  frequen¬ 
cies  below  0.3  rad/sec  are  due  to  the 
phugold  dynamics. 

The  coherence  function  shown  in 
Fig.  9  la  strong  over  the  frequency 
range  of  Interest  (0.2-6 .0  rad/sec), 
as  shown.  For  Input  frequencies  above 
7.0  rad/sec,  the  coherence  function 
becomes  erratic  and  the  transfer- 
function  Identification  Is  less  accu¬ 
rate.  For  low-frequency  Inputs  (less 
than  0.3  rad/sec),  the  pitch-rate 
response  decreases,  even  for  the 
nearly  constant  Input  amplitude  owing 
to  the  effect  of  the  phugold  dynamics. 

This  results  In  a  decrease  of  Informa¬ 
tion  transfer  and  an  associated  drop 
in  the  coherence  function.  This 
coherence-function  roll-off  is  also 
attributable  to  atmosphere  turbulence 
effects,  which  become  more  Important 
at  low  frequency  where  the  turbulence 
spectrum  peaks  (Ref.  [7]). 

The  Bode  plot  of  Figs.  8(a)  and  8(b)  Is  the  basic  format  for  presenting  classical 
frequency-domain  results.  It  provides  useful  information  on  dynamic  characteristics  and  on 
handling  qualities  Implications.  To  the  extent  that  the  linear-describing  function  Is  an  accu¬ 
rate  representation  of  the  Input-to-output  process,  the  Bode  plot  is  an  exact  description  of  the 
vehicle  dynamics  and  contains  no  Implicit  assumptions  on  model  order  or  structure.  This  Is  the 
key  advantage  of  frequency-domain  Identification  over  the  more  conventional  time-domain 
approaches . 

Recent  military  specifications  for  plloted-handling  qualities  (Ref.  [4])  are  based  on 
two  key  frequency-domain  parameters:  bandwidth  (ug^f)  and  effective  time  delay  (tp)  .  For  a  par¬ 
ticular  transfer  function,  the  bandwidth  Is  defined  as  that  frequency  where  the  phase  margin  Is 
45  deg,  or  the  gain  margin  Is  6  dB,  whichever  Is  lowest.  As  discussed  in  Ref.  [4],  the  band¬ 
width  Is  a  measure  of  the  speed  of  response;  a  high  bandwidth  reflects  quick  response  and  accu¬ 
rate  tracking  capability;  a  low  bandwidth  suggests  sluggish  response  and  pilot-induced  oscilla¬ 
tion  tendencies.  The  effective  time  delay  Is  a  measure  of  the  slope  of  the  phase  curve  for 
frequencies  near  the  bandwidth  value.  A  large  effective  time  delay  (e.g.,  greater  than 
100  msec)  results  In  a  significant  degradation  of  piloted  handling  qualities. 

For  longitudinal  control,  a  key  transfer  function  is  pitch  attitude  response  to  elevator, 
e/6g,  which  Is  easily  derived  from  the  pitch-rate  transfer  function  as: 

9(8)  .  q(») 

«^(s)  86^(8) 


Figure  5.  Pitch-rate  response  (q)  during  longitudinal 
frequency  sweeps . 


Figure  6.  Pitch-rate  output  autospectrum  (Gq^)  . 


Figure  7.  Magnitude  of  cross-spectrum  between  ele¬ 
vator  surface  Inputs  and  pitch-rate  response  (|G{  q|). 


(4) 


Figure  8.  Pitch-rate  response  to  elevator  (q/-5g) 
In  cruise,  (a)  Transfer  function  magnitude. 

(b)  Transfer  function  phase. 


Figure  9.  Coherence  function  for  pitch-rate 
response  Identification  (Yj^q) • 


-  FLIGHT  DATA 

- TRANSFER  FUNCTION  MODEL 


Figure  10.  Vertical-acceleration  response  to 
elevator  (a^/d.)  in  cruise,  (a)  Transfer-function 
magnitude,  (b)  Transfer-function  phase. 


where  s  is  the  Laplace  operator. 

The  Bode  plot  for  pitch-attitude 
response  may  be  obtained  from  the 
pitch-rate  response  by  rotating  the 
magnitude  curve  about  the  u  -  1  rad/sec 
point  by  slope  of  20  dB/decade.  The 
phase  curve  is  shifted  by  a  constant 
value  of  -90  deg.  Referring  to  Fig.  8, 
the  critical  bandwidth  criterion  is 
phase  margin  which  Is  45  deg  at  a  fre¬ 
quency  of  "2.5  rad/sec.  The 

effective  time  delay  is  negligible  since 
the  phase  curve  is  nearly  flat  where  it 
is  near  180  deg.  These  bandwidth  and 
time  delay  values  are  well  within  the 
Ref.  [4]  Level  II,  Category  A,  specifi¬ 
cations  (adequate  handling  qualities  for 
high-precision  tracking  tasks) ,  and  cer¬ 
tainly  reflect  satisfactory  character¬ 
istics  for  SCAB  failure  conditions. 

The  transfer  function  for  novmaL 
acceleration  ree-ponee  to  elevator, 
a^/dg,  is  shown  in  Figs.  10(a)  and  10(b). 
The  response  is  dominated  by  the  classi¬ 
cal  second-order  short-period  mode  over 
the  majority  of  the  frequency  range. 

The  magnitude  curve  is  flat  at  mid- 
frequency  indicating  a  constant  normal 
acceleration  response  to  a  step  elevator 
input,  with  a  roll-off  in  response  for 
frequencies  beyond  the  short  period  mode. 
The  fall-off  in  normal  acceleration 
response  for  frequencies  below 
0.3  rad/sec  is  probably  due  to  the 
dominance  of  the  phugold  dynamics.  The 
phase  curve  exhibits  the  classical  second 
order  response;  l.e.,  zero  phase  lag  at 
low  frequency,  180  deg  of  phase  lag  at 
high  frequency,  and  90  deg  of  phase 
lag  at  the  second  order  mode 
(u  -  2.0  rad/sec).  This  is  consistent 
with  the  previous  pitch-rate  results. 

The  coherence  function  f'~r  the  normal 
acceleration  response  shown  in  Fig.  11 
is  strong  over  the  frequency  range  of 
0.3  to  10.0  rad/sec,  with  the  fall-off 
at  low  frequency,  again  owing  to  the 
dominance  of  phugoid  dynamics  and  tur¬ 
bulence  effects.  As  before,  this  sug¬ 
gests  excellent  identification  over  the 
entire  frequency  range  of  interest. 

Longitudinal  Transfer  FYinction  Fitting 

Analytical  transfer  function 
forms  are  selected  for  each  degree  of 
freedom  based  on  configuration  and  flight 
condition  factors.  In  the  hover  flight 
condition,  vehicle  dynamics  are  dominated 
by  the  hovering  cubic,  and  decoupled 
heave  and  yaw  modes.  In  wing-borne 
flight,  the  conventional  longitudinal 
and  lateral  quartlc  equations  dominate. 
Thus,  transfer-function  models  which  are 
appropriate  for  the  hover  flight  condi¬ 
tion  are  not  necessarily  applicable  to 
forward  flight  conditions.  Obviously,  if 
a  model  of  high  enough  order  is  selected, 
the  parameters  can  be  adjusted  to  accom¬ 
modate  each  flight  condition.  However, 
such  transfer  function  models  no  longer 
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retain  the  physical  significance  of  the 
classical  lower-order  parameters. 

Also,  higher-order  transfer-function 
models  tend  to  be  strongly  tuned  to  the 
specific  Inputs  which  are  used  In  the 
Identification  procedure  (e.g.,  a  fre¬ 
quency  sweep)  and  are  often  very  poor 
predictors  of  other  test  Inputs  (e.g., 
step  Inputs) .  Therefore,  higher-order 
models  are  not  desirable.  In  the 
approach  taken  In  Ref.  [3],  the 
minimum-order  transfer-function  models 
'which  can  satisfactorily  fit  the'  Fre¬ 
quency  responses  are  used  with  the 
upper  limit  taken  as  the  physical  order 
of  the  system. 


-  -  FREQUENCY,  rad/*Bc- 

Flgure  11.  Coherence  function  for  vertical- 
acceleration  response  Identification  ’ 


Examination  of  the  longitudinal  frequency  responses  of  Figs.  8(a)  and  8(b)  shows  that  the 
longitudinal  dynamics  of  this  flight  condition  are  dominated  by  the  short  period  mode 
(ugp  -  2  rad/sec).  Therefore,  we  adopt  the  classical  pitch  rate  and  normal  acceleration 
responses  to  elevator*: 

q(s)  _  ^ 


where  q(s)/6_(s)  Is  thi>  Laplace-transformed  pitch  rate  response  to  elevator  surface  deflection, 
deg/sec/deg-elevator ;  Is  the  elevator  pitch  sensitivity;  1/Tg  Is  the  first-order  numera¬ 


tor  Inverse  time  constant ; 


are  the  equivalent  short-period  mode  damping  and  natural 


- - - ..gp,  -,gp  - - - 

frequency,  respectively;  and  T0  Is  the  effective  tine  delay;  and 


f5sp-  “spl 

where  a^(,a) Is  the  Laplace-transformed  vertical  acceleration  response  (positive  downward) 
to  elevator  surface  deflection,  g/deg-elevator ;  is  the  elevator  vertical  sensitivity;  and 

the  denominator  parameters  are  identical  to  those  of  Eq.  (5).  The  effective  time  delay  for  the 
vertical  acceleration  response  is 

We  Ignore  the  low-frequency  phugold  dynamics  since,  as  seen  from  Figs.  8  and  1C,  these 
are  important  only  for  the  very  lowest  frequency  Inputs. 

The  longitudinal  transfer  functions  of  Eqs.  (5)  and  (6)  have  the  same  denominator  factors, 
so  simultaneous  fitting  of  these  responses,  while  imposing  the  restriction  of  commonality  of 
denominator  parameters,  makes  good  physical  sense.  Also,  the  analyses  of  Refs.  [3]  and  [8], 
suggest  that  such  simultaneous  fitting  techniques  are  needed  to  ensure  unique  and  physically 
meaningful  values  for  the  numerator  factors.  For  the  longitudinal  transfer  functions,  this 
avoids  the  otherwise  so-called  "fixed  or  free"  Lg  problem  (Ref.  [4]). 

The  transfer-function  parameters  of  Eqs.  (5)  and  (6)  are  Iteratively  varied  in  order  to 
obtain  a  best  least-squares  fit  between  the  equations  and  the  frequency  responses  of  Figs.  8 
and  10  over  the  selected  frequency  range.  This  procedure  is  completed  using  the  computer  pro¬ 
gram  UNF^  developed  by  Hodgklnson  et  al.  (Ref.  [9]).  For  the  present  case,  the  selected  range 
for  simultaneous  fitting  of  the  pitch  rate  and  normal  acceleration  responses  is  0.3  to  7  rad/sec. 
In  this  frequency  range,  the  dynamics  are  clearly  dominated  by  the  short-period  mode  and  the 
coherence  Is  strong  for  both  transfer  functions.  Once  the  short-period  damping  and  natural  fre¬ 
quency  (Csp,  o>gp)  are  obtained  for  Che  simultaneous  fit,  the  high-frequency  gain  for  the  normal 
acceleration  response  (Z^^)  Is  varied  holding  the  daisping  and  frequency  constant.  This  opti¬ 
mizes  the  a2  fit  over  Che  frequency  range  0.3  to  10.0  rad/sec  since  Che  coherence  function  for 
this  measurement  remains  strong  out  to  higher  frequencies.  The  transfer  function  parameters  for 
Che  pitch  race  and  normal  acceleration  responses  are  finally  obtained  as: 


*Shorthand  notation:  [;,  u]  Implies  s’  ■¥  ZCus  -f  u’,  (  •  damping  ratio,  u  -  undamped 
natural  frequency  (rad/sec);  and  (1/T)  implies  s  +  (1/T),  rad/sec. 


(7) 


Mr  =  -1  .in  deg/sec*/deg-elevator^ 
“e 

-  -32.144  deg/8ec^/ln.-6Lou 

l/Tft  =  1.035  rad/sec 
”2 

?sp  =  0.554 

oj__  =  2.179  rad/sec 

sp 

Tg  “  0.016  sec  __ 


-Zfe-- 


4 


sp 


u 


sp 


T 


‘Z 


1.597  gnleg-elevacot 
6.644  g/in.-6L0{j 

0.554 

2.179  rad/sec 
0.018  sec 


(8) 


A  comparison  of  these  lower-order  models  with  Che  flight  test  results  Is  presented  In  Figs.  8 
and  10.  The  pitch-response  m'<^l  matches  Che  data  very  well  over  the  selected  fitting  range 
(0.3-10.0  rad/sec  for  q/6p  ,  wiilch  shows  that  the  short  period  approximation  adequately  repre¬ 
sents  the  high-  and  mld-f.r  uency  dynamics.  For  low-input  frequencies  (u  <  0.3  rad/sec),  the 
phugold  mode  causes  a  drop  In  the  pitch  response,  a  characteristic  which  is  not  "captured"  by 
Che  short  period  model.  The  match  between  the  normal  acceleration  transfer  function  and  the 
flight  data  Is  not  nearly  as  good  In  Che  fitting  range,  and  shows  a  noticeable  discrepancy  In 
the  magnitude  response  for  frequencies  of  greater  than  2  rad/sec.  Possible  sources  of  the  mid- 
and  high-frequency  discrepancy  are  the  rotor  rpm  dynamics,  which  are  known  to  be  important  for 
the  forward  flight  conditions  (Ref.  [10])  and  rotor  inflow  dynamics.  Neither  of  these  sources 
of  added  dynamics  Is  directly  accounted  for  by  the  classical  short-period  approximation,  except 
to  the  extent  that  the  short  period  parameters  are  adjusted  in  order  to  best  fit  the  flight 
data.  As  in  the  pitch  response,  the  low  frequency  mismatch  is  due  to  the  omission  of  the  phu¬ 
gold  mode  from  the  transfer-function  model. 

Guidelines  are  given  in  the  Mil  Handbook  (Ref.  [4])  for  the  maximum  allowable  mismatch  in 
lower-order  system  modeling.  The  discrepancies  shown  In  the  normal  acceleration  responses 
(Figs.  10(a)  and  10(b))  are  well  within  these  allowable  guidelines,  so  the  identified  parameters 
should  be  useful  Indicators  of  plloted-handllng  qualities.  The  capability  to  directly  evaluate 
the  effect  of  model  order  on  response  fitting  Is  clearly  seen  in  this  example,  again  represent¬ 
ing  a  unique  advantage  of  the  frequency  domain  approach.  Further  analyses  could  include  the 
rotor  rpm  and  phugold  degrees-of-freedom  in  the  transfer-function  models,  and  could  evaluate 
their  Importance  with  respect  to  vehicle  dynamics  and  handling  qualities. 

Lower  order  system  parameters  are  useful  in  characterizing  the  vehicle  dynamics  for  com¬ 
parison  with  the  dynamics  of  other  aircraft.  The  identified  short  period  damping  and  frequency 
(Cgp  =  0.554  and  Ugp  =  2.18  rad/sec)  are  consistent  with  preliminary  observations  which  are 
based  on  the  raw  Bode  plot  Information.  The  relatively  high  degree  of  damping  is  a  reflection 
of  the  small  response  peak  of  the  pitch-rate  transfer  function.  The  small  effective  time  delays 
(In  the  pitch  attitude  and  normal  acceleration  transfer  fits)  Indicate  that  the  high  frequency 
flexibility  and  servo-lag  effects  are  not  Important  for  this  flight  condition,  as  noted  earlier. 
The  Mil  Handbook  requirements  for  lower-order  equivalent  system  pitch  response  are  given  In 
terms  of  the  parameters  UgpTej,  4sp,  and  tg.  The  values  of  these  parameters  given  In  Eqs.  (7) 
and  (8)  are  well  within  the  Level  I,  Category  A,  handling  qualities  requirements,  confirming 
that  the  longitudinal  characteristics  for  SCAS  off  conditions  are  adequate. 

Time-domain  verification  of  the  longitudinal  transfer  function  models  will  be  addressed 
after  lateral  response  models  are  presented. 

Lateral  Dynamics 

The  open-loop  transfer  function  for  roll-rate  reeponee  to  aHerona  p/6g,  is  shown  in 
Figs.  12(a)  and  12(b).  This  Is  a  classical  first-order  system  in  which  aileron  Inputs  produce 
a  constant  roll  acceleration  at  high  frequency,  and  a  constant  roll  rate  at  low  frequency.  The 
corner  frequency  Is  about  1  rad/sec,  with  an  associated  phase  lag  of  roughly  -45  deg,  as 
expected.  The  dominant  time  constant  is  thus  roughly  1  sec.  Implying  about  2  to  3  sec  to  reach 
a  steady-state  roll  rate.  The  coherence  function  shown  In  Fig.  13  is  very  strong  over  the  fre¬ 
quency  range  of  0.1  to  9.0  rad/sec.  Indicating  an  excellent  identification  of  the  roll-rate 
dynamics. 
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Specifications  for  the  lateral 
handling  qualities  of  military  air¬ 
craft  are  given  In  the  Mil  Handbook 
(Ref.  [4]),  based  on  a  lower-order 
equivalent-system  model  for  the  roll- 
attitude  dynamics  (♦/Sg).  When  roll/ 
yaw  coupling  is  minor,  the  specifica¬ 
tion  on  the  equivalent  system  roll 
mode  time  constant  (Tj)  can  be  Inter¬ 
preted  as  the  reciprocal  of  the  roll 
response  bandwidth  (ugy^) .  Referring 

again  to  Fig.  12(b),  the  bandwidth  for 
'fbri:' attitude  response  to  ailerons 'Is 
“BW^  =0.9  rad/sec,  yielding  a  domi¬ 
nant  time  constant  of  about  1.1  sec. 
This  satisfies  the  Ref.  [4]  specifica¬ 
tion  (Tp  <1.4  sec)  for  Level  I,  Cate¬ 
gory  A,  handling  qualities. 

Aerodynamic  sideslip  measure¬ 
ments  are  obtained  from  a  sideslip 
Indicator  located  roughly  18  ft  ahead 
of  the  aircraft  c.g.  The  sideslip  at 
the  c.g.,  Bpg,  Is  calculated  by  cor¬ 
recting  the  measured  signal  for  posi¬ 
tion  error  based  on  yaw  rate  and  air¬ 
speed.  No  corrections  are  made  for 
sensor  dynamics  because  these  are  not 
felt  to  be  significant  within  the 
Identification  bandwidth. 

The  transfer  function  for  side¬ 
slip  response  (at  the  c.g.)  to  rudder 
inputs  Bcg/^ri  fs  shown  in  Figs.  14(a) 
and  14(b).  The  response  is  character¬ 
ized  by  a  lightly  damped  second-order 
mode  with  a  frequency  of  about 
1.6  rad/sec.  High-frequency  rudder 
Inputs  yield  a  constant  sideslip  accel¬ 
eration,  while  low-frequency  rudder 
Inputs  yield  a  constant  sideslip  angle 
corresponding  t.)  a  steady-state  yaw 
rate.  The  coherence  function  for  this 
case,  shown  in  Fig.  15  is  strong  over 
the  frequency  range  0.1  to  5.0  rad/sec, 
falling  off  sharply  in  the  high  fre¬ 
quency  range. 

The  rapid  and  relatively  low- 
frequency  decline  of  the  coherence 
function  of  Fig.  15  is  unlike  the  pre¬ 
vious  transfer  functions  considered 
for  the  cruise  flight  condition.  The 
reason  for  this  difference  is  that 
aerodynamic  sideslip  Is  a  state  vari¬ 
able  which  Is  one  derivative  lower 
than  the  other  rate  variables  consid¬ 
ered  previously.  The  response  of 
parameters  such  as  roll  rate  and 
pitch  rate  falls  off  at  high  fre¬ 
quency  with  a  slope  of  20  dB/decade. 
However,  the  aerodynamic  sideslip 
response  as  seen  In  Fig.  14(a)  falls 
off  at  the  much  faster  rate  of 
40  dB/decade.  Although  the  Input 
autospectrum  Is  roughly  constant 
over  the  entire  plotted  frequency 
range,  this  rapidly  falling  response 
causes  a  marked  reduction  of  the 
output  autospectrum  and  cross¬ 
spectrum  and  an  associated  decline 
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Figure  12.  Roll-rate  response  to  aileron  (p/fig)  in 
cruise,  (a)  Transfer-function  magnitude. 

(b)  Transfer-function  phase. 
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Figure  13.  Coherence  function  for  roll-rate  response 
identification  (y^  p) . 
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Figure  14.  Sideslip  response  to  rudder  (B(,g/6j.)  In 
cruise,  (a)  Transfer-function  magnitude. 

(b)  Transfer-function  phase. 


Figure  IS.  Coherence  function  for  sideslip  response 
Identification  • 


In  the  coherence  function.  This  enpha-  i 

sizes  the  advantages  of  using  the 
higher  state  derivatives  In  Identify¬ 
ing  the  high-frequency  dynamics  and 
suggests  the  possible  advantages  of 
using  attitude  signals  to  Identify  the 
low-frequency  dynamics.  However,  since 
the  yaw  mode  of  Interest  has  a  natural 
frequency  of  roughly  1.6  rad/sec.  the 
falling  coherence  for  frequencies 
greater  than  5  rad/sec  Is  not  a  severe 
limitation. 

As  In  the  roll  axis,  handling 
qualities  specifications  for  the  yaw 
axis  response  are  given  In  terms  of 
equivalent  system  models.  These  are 
derived  In  the  following  section. 


Lateral  Transfer  Function  Fitting 

Equivalent  system  fitting  using  decoupled  and  coupled  lateral/directional  models  has 
been  considered  in  detail  by  Blschoff  (Ref.  [8]).  Uith  the  deooupled  model,  the  responses  are 
fit  independently,  using  a  first-order  roll  rate  and  a  second-order  sideslip  transfer  function 
The  coupled  model  approach  Is  based  on  simultaneous  fitting  of  the  roll  and  sideslip  responses 
to  obtain  the  following  fourth-order  transfer  functions: 


P(s)  * 

«a(s)  (l/Tg)(l/Tp[Cd.  “d^ 


where  p(8)/dg(s)  Is  the  Laplace-transformed  roll  rate  response  to  aileron  surface  deflection, 
deg/sec/deg-aileron;  Lj^  Is  the  aileron  roll  sensitivity;  are  the  second-order  numera¬ 

tor  damping  and  natural  frequency,  respectively;  I/T3,  l/T^  are  the  equivalent  spiral  and  roll 
subsidence  modes;  are  the  equivalent  Dutch  roll  mode  damping  and  natural  frequency, 

respectively;  and  Is  the  effective  time  delay. 

6„(«) 

«,(s)  “  (l/T^)(l/VUd'  “dl  ^  ' 


where  8cg(s)/£f (s)  Is  the  Laplace-transformed  sideslip  response  to  rudder  surface  deflection, 
deg/deg-elevator ;  Is  the  rudder  sideslip  sensitivity;  1/Tgj,  l/Tg^,  l/Tgj  are  the  first- 

order  numerator  Inverse  time-constants;  and  the  denominator  parameters  are  identical  to  those 
of  Eq.  (9).  The  effective  time  delay  for  the  sideslip  response  Is  Tg.  The  simultaneous  fit¬ 
ting  approach  is  consistent  with  that  used  In  the  longitudinal  axis  and  allows  the  Identifica¬ 
tion  of  the  u./uj  "coupling  effect"  In  the  roll  response  (Eq.  (9),  idiich  is  iaiportant  In  the 
handling  qualities  assessment. 


The  coherence  function  results  of  Figs.  13  and  15  show  a  satisfactory  Identification  of 
both  roll  and  sideslip  responses  in  the  frequency  range  of  0.1  to  5.0  rad/sec.  Based  on  this 
fitting  range,  the  computer  program  LATFIT  (Ref.  [11])  was  used  to  obtain  the  paraoieters  of 
Eqs.  (9)  and  (10),  sianiltaneously .  Since  the  coherence  for  the  roll  reeponee  Is  satisfactory 
over  the  entire  range  of  0.1  to  10.0  rad/sec,  this  degree  of  freedom  was  refit  alone,  holding 
the  denominator  factors  constant  at  the  values  obtained  from  the  simultaneous  solution.  This 
procedure  Is  analogous  to  that  used  In  the  longitudinal  case  and  optimizes  the  values  of  the 
high-frequency  gain  and  numerator  parameters  In  the  roll  transfer  function.  The  final  results 
for  the  parameters  of  Eqs.  (9)  and  (10)  are: 


l/Tg 

1/T^ 

5d 


4.486  deg/sec^/deg-alleron  "N 
17.630  deg/8ec*/ln.-6u^.j 


0.313 

1.887  rad/sec 
0.063  rad/sec 
1.090  rad/sec 
0.248 


(11) 
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The  lateral/direct lonal  cranafer  functian  flta  are  plotted  aa  dotted  llnea  for  conpari- 
son  with  the  flight  data  of  Figs.  12  and  14.  The  watching  of  the  roll-response  siagnltude 
(Fig.  12)  Is  excellent  over  the  entire  frequency  range.  Only  a  alight  anoisaly  In  the  phase 
watching  Is  apparent  for  low  frequency  Inputs.  The  agreewent  between  the  sideslip  fit  and  the 
flight  data  la  alao  excellent  (Fig.  14). 

As  a  result  of  the  slaniltaneous  fitting  of  the  roll  and  yaw  responses >  the  transfer 
function  paraweters  of  Eqs.  (9)  and  (10)  have  clearly  retained  their  physical  significance. 

The  Inverse  roll  mode  time  constant.  1/T,  •  1.090  rad/sec  la  roughly  equal  to  the  roll  responae 
bandwidth  aa  expected;  both  are  well  within  the  level  I  roll  response  requlrewents,  as  wentloned 
earlier.  The  equivalent  Dutch  roll  node  la  lightly  dawped  with  a  natural  frequency  roughly 
corresponding  to  the  peak  in  the  sideslip  Bode  wagnltude  plot  (Fig.  14(a)).  The  swall  effective 
tiwe  delays  for  the  roll  and  sideslip  responses  reflect  negligible  high-frequency  flexibility 
and  servo-lag  effects,  and  are  consistent  with  the  previous  longitudinal  results.  The  swall 
tiwe  delay  In  the  sideslip  response  also  supports  the  owission  of  corrections  for  sideslip 
sensor  dynawlcs. 

The  ratio  of  the  natural  frequency  of  the  nunerator  complex  zero  (u^)  to  the  denowlnator 
cowplex  pole  (uj)  Is  one  measure  of  roll/yaw  coupling  in  response  to  aileron  Inputs.  When 
•  uj,  the  numerator  and  denominator  quadratic  factors  roughly  cancel,  and  the  resulting 
decoupled  roll  response  is  characterised  entirely  by  the  roll-mode  time  constant.  This  case 
leads  to  the  best  handling  qualities  for  a  nominal  value  of  the  roll-mode  tiwe  constant.  As 
the  roll/yaw  coupling  Increases,  the  numerator  and  denominator  quadratic  factors  of  Eq.  (9)  no 
longer  cancel  and  an  undesirable  oscillatory  component  of  roll  rate  la  generated.  Referring  to 
Eq.  (11),  the  near  unity  value  -  1.19  suggests  no  such  concern  for  roll/yaw  coupling. 

A  iBore  direct  measure  of  the  roll/yaw  coupling  Is  obtained  from  evaluating  the  ratio  of 
roll  attitude  to  sideslip  evaluated  at  the  Dutch  roll  frequency  |4/6|j.  Analyses  show 
(Ref.  [4])  that  there  Is  a  direct  correlation  between  the  Id/Blj  ratio  and  the  ratio  of  oscil¬ 
latory  roll  rate  to  ateady  roll  rate  for  step  aileron  Inputs.  A  small  value  of  |d/B|d 
less  than  1.4)  suggests  a  small  oscillatory  roll-rate  component  and  desirable  handling  quali¬ 
ties.  A  large  value  causes  roll  response  overshoot  and  related  poor  tracking  characteristics. 

A  cross  correlation  of  roll  rate  and  sideslip  responses  yields  |4/B|d  **  1.27,  idilch  Implies  no 
such  tracking  deficiencies.  The  preceding  analyses  have  shown  that  the  coll  response  Is  largely 
decoupled  and  first  order  In  nature.  Therefore,  the  slaiple,  first-order  roll  rate  model  could 
be  adopted  for  future  studies. 


-  1.581  rad/sec 
”  0.045  sec  7 

-0.051  deg/sec/deg-rudder' 
-0.408  deg/sec/ln.-Bpgj 

0.086  rad/sec 
0.818  rad/ sec 
47.946  rad/sec 
0.063  rad/sec 
1.090  rad/sec 
0.248 

1.381  rad/sec 
0.026  sec 


(11) 

(cont.) 


(12) 


4.  Tlsie  Domain  Verification  of  Transfer  Function  Hodela 


The  transfer  function  is  the  minimum  realisation  description  of  a  linear  Input-output 
process.  The  transfer-function  models  developed  above  can  be  used  to  generate  all  other  fre¬ 
quency  and  time  domain  Information.  Therefore,  a  set  of  verified  transfer  function  aK>dels  la  a 
useful  output  format  for  the  linear  Identification  procedure.  However,  the  adequacy  of  these 
linearised  transfer-function  models  for  predicting  the  tlme-doMln  dynawlcs  of  nonlinear  systems 
la  at  question.  For  such  nonlinear  systems,  the  transfer  function  Is  actually  a  describing 
function  which  Is  strictly  valid  only  for  the  Input  amplitudes  which  were  used  In  the  flight 
test  experiment.  Also,  since  Identification  procedures  tend  to  be  tuned  to  the  type  of  test 
Inputs  used  (e.g.,  a  sine-sweep),  their  accuracy  In  predicting  responses  to  other  classes  of 


Inputs  (e.g..  steps,  doublets)  is  often  uncertain.  Lastly,  although  the  aismatches  observed  in 
the  previous  section  between  the  transfer-function  fits  and  the  frequency-response  flight  data 
never  exceeded  the  allowable  handling  quality  specifications,  the  laportance  of  these  discrep¬ 
ancies  in  predicting  tiae  response  behavior  is  still  at  issue.  These  questions  are  addressed 
in  the  following  section. 

The  basic  verification  approach  is  to  conpare  the  step  response  of  the  aircraft  with  the 
response  of  the  transfer-function  model  driven  by  the  same  recorded  control  inputs.  The  step 
responses  were  cosq>leted  for  all  flight  conditions  with  the  stabilization  system  disengaged  in 
order  to  show  that  the  extracted  transfer-function  models  are  valid  for  the  bare-airfraae  con¬ 
figuration.  The  transfer  functions  for  the  hover  flight  condition  and  those  extracted  in  the 
present  analysis  for  the  cruise  flight  condition  are  sunmarlzed  in  Table  1.  For  illustrative 
consistency,  all  transfer  functions  are  referenced  to  the  three  surface  deflections:  elevator, 
aileron,  and  rudder.  The  verification  results  for  the  hover  flight  condition  are  presented 
first . 


TABLE  1.-  SUMMARY  OF  IDENTIFIED  TRANSFER-FUNCTION  MODELS^  FOR  THE 
HOVER  AND  CRUISE  FLIGHT  CONDITIONS 


Hover  (V  =  0) 


Cruise  (V  ■  170  knots,  indicated) 


Longitudinal  dynamics 


-0.108  e"*'*"*® 


(0.115) 


;  g/in.—  power  lever 


S  1.597  e~*-*^*^ 
“  {0.554,  2.179) 


:  g/deg- elevator 


-2.607s(-0.075)(0.686)e 
«  "  ■ 


'0.0I|4S 


(0,115)(1.597)(-0.978,  0.393)  ’ 

deg/sec/deg  —  elevator 


-1. 

«e 


-7.727(1.035)e 


-0 .01(8 


(0.554,  2.179) 


deg/sec/deg—  elevator 


Lateral  dynamics 


JL 

6. 


r 

6_ 


r 

6. 


-0 .022S 


3.959a(-0.092)(0.570)e _ 

(0.108) (1.429) [-0.581,  0.427]  ‘ 

deg/sec/deg  —  aileron 


0.708  e-***^*° 
(0.108) 


deg/sec/deg—  rudder 


«r 


4.486s[0.313.  1.887]e~* 

(0.063) (1.090) [0.248,  1.581)  ' 


deg/sec/deg  —  aileron 


-0s02»8 


0.398 (-0.310) [0.762,  0.524]e~° •***^  . 
(0.108) (1.429) [-0.581,  0.427) 

deg/sec/deg—  aileron 


-0.051(0.086) (0.818) (47.946)e 

(0.063)(1.090)[0.248,  1.581)  ‘ 

deg/deg  —  rudder 


‘^Shorthand  notation:  [c,  id]  implies  s^  24(ds  -i-  ^  •  damping  ratio,  id 

frequency  (rad/sec);  and  (1/T)  implies  s  +  (1/T),  rad/sec. 


undamped  natural 


Hover  Flight  Condition 

A  co-plot  of  the  transfer-function  model  for  yaw  rate  response  to  rudder  r/6-  (Table  1) 
with  the  frequency  response  data  of  Ref.  [3]  is  presented  in  Fig.  16.  The  transfer-function 
model  matches  the  flight  data  very  well  at  high  frequency  in  both  magnitude  and  phase  but  not  as 
well  at  low  frequency.  Even  so,  the  fit  meets  the  mismatch  criteria  of  Ref.  (4)  for  handling 
qualities  assessment.  These  criteria  emphasize  the  piloting  crossover  range  of  1-3  rad/sec,  in 
lAich  the  fit  is  excellent  for  this  case. 


When  the  transfer-function  model  is  driven  with  one  of  the  swept  sine-woue  inputs  used  in 
the  identification  procedure,  the  time-domain  ramifications  of  model  mismatch  are  clearly 
exposed.  The  transfer-function  model  and  flight-data  responses  for  a  typical  frequency  sweep 
are  compared  in  Fig.  17.  The  input  amplitude  shown  in  Fig.  17(a)  is  roughly  i20Z  of  the  maximum 
pedal  travel  producing  a  yaw  response  amplitude  of  roughly  ±20  deg/sec,  as  seen  in  Fig.  17(b). 
The  flat  spot  in  the  yaw  rate  signal  at  50  sec  is  due  to  yaw-rate  sensor  limiting.  In  general, 
the  Mtch  between  the  transfer-function  response  and  the  flight  data  is  excellent.  The  phasing 
and  form  of  the  model  response  matches  the  flight  data  very  well,  with  the  only  noticeable 
anomaly  being  an  underestimation  of  the  peak  amplitudes  for  low-frequency  Inputs.  This  is  a 
result  of  the  fit  mismatch  at  low  frequency  (Fig.  16).  Thus,  the  transfer  function  model  accu¬ 
rately  captures  the  essence  of  the  yaw-rate  response  over  the  frequency  range  of  Interest,  and 
can  be  considered  adequate  for  sine-wave  type  Inputs. 
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The  flight  data  for  yaw  rate 
response  to  a  step  pedal  input  is  shown 
in  Figs.  18(a)  and  18(b).  The  control 
deflection  represents  roughly  40Z  of 
the  maximum  control  power,  a  fairly 
large  amplitude  input.  The  maximum 
yaw-rate  response  is  about  35  degVsec. 
«dilch  is  75Z  greater  than  the  peak 
values  obtained  in  the  frequency 
sweeps.  The  transfer-function  model 
response  co-plotted  in  the  dashed 
curve  of  Fig.  18(b)  compares  extremely 
well  with  the  flight  data  over  the 
entire  run—  even  for  this  fairly  large 
input  <UBplltude  case.  Analyses  on  the 
vertical  response  to  collective 

showed  the  same  excellent 
response  correlation.  This  Is 
expected  since  the  response  dynamics 
for  the  heave  and  yaw  axes  are  nearly 
the  same  (Table  1) . 

The  step-response  run  lengths  of 
roughly  10  sec  emphasize  the  high-  and 
mid-frequency  ranges,  allowing  only 
enough  time  for  about  one  time  constant 
(Ty  ”  9.3  sec).  Thus,  the  low- 
frequency  mismatches  never  really  have 
enough  time  to  build  up.  This  is 
typical  of  step-response  data  and 
Illustrates  the  advantages  of  a  sym¬ 
metrical  alne-sweep  input.  Clearly, 
however,  the  response  In  the  first  few 
seconds  la  the  key  concern  of  the 
pilot.  This  Is  reflected  In  the  mis¬ 
match  criteria  being  most  stringent  In 
the  mid-frequency  range. 

So  far,  the  time  response 
results  show  that  the  simple  transfer- 
function  models  match  the  step-response 
characteristics  very  well.  Also,  the 
use  of  a  very  low-order  description  has 
avoided  any  problems  In  over-tuned 
modeling.  In  general,  however,  good 
time  history  matching  is  highly  depen¬ 
dent  on  the  level  of  spurious  Inputs 
to  the  model  due  to:  (1)  Improper 
setting  of  Initial  conditions, 

(2)  control  Inputs  required  to  sup¬ 
press  turbulence,  and  (3)  unmodeled 
control  coupling.  When  the  model 
dynamics  are  stable,  the  effects  of 
these  spurious  Inputs  die  out  with 
time,  and  response  matching  presents 
no  particular  problems.  Such  cases 
are  the  easiest  for  time-domain  Iden¬ 
tification  as  well. 


Figure  16.  Yaw-rate  response  to  rudder  (r/S^)  in 
hover,  (a)  Transfer-function  magnitude. 

(b)  Transfer-function  phase. 


-  FLIGHT  DATA 

- TRANSFER  FUNCTION  MODEL 

5.-kw\/^ 


-20- 


<■) 
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Figure  17.  Comparison  of  aircraft  and  transfer- 
function  model  responses  to  a  pedal  frequency-sweep 
in  hover,  (a)  Rudder  surface  deflection,  (b)  Yaw- 
rate  response. 


When  the  dynamics  of  the  transfer-function  response  are  unstable,  the  transients  due  to 
these  spurious  Inputs  grow  with  time.  Therefore,  small  errors  in  the  choice  of  Initial  condi¬ 
tions,  or  small  Inputs  required  to  suppress  the  response  to  turbulence,  can  create  gross  differ¬ 
ences  between  the  transfer- function  response  and  the  flight  data.  Also,  for  unstable  vehicles, 
acquisition  and  maintenance  of  a  steady  trim  condition  Is  often  not  possible,  further  coaq>llcat- 
Ing  the  proper  initialization  of  the  transfer  function  calculations.  As  the  level  of  Instabil¬ 
ity  of  the  transfer-function  dynamics  Increases,  so  do  the  problems  of  time-history  matching. 

Mot  surprisingly,  these  are  also  the  worst  cases  for  time-domain  Identification. 

The  roll-rate-to-alleron  transfer  function  (Table  1)  is  dominated  by  an  unstable  roll 
oscillation  having  a  time-to-double  of  about  2.8  sec.  This  reflects  the  speed  at  Milch  the 
errors  and  Inexact  control  inputs  and  Initial  conditions  will  propagate.  A  good  step  response 
In  the  roll  axis  was  achieved  starting  from  a  fairly  steady  Initial  trim  condition  as  shown  In 
Figs.  19(a)  and  19(b).  The  assumed  Initial  condition  Is  determined  from  the  average  of  the 
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•lop«>  rise  tlae,  and  percent  overshoot 
are  all  well  predicted;  therefore,  the 
transfer-function  andel  is  a  useful 
characterization  of  the  acceleration 
response.  This  cosparison  further  indi¬ 
cates  the  need  to  investigate  sources 
of  unnodeled  dynamics  (e.g..  engine 
response)  in  the  vertical  axis. 

The  aircraft  roll  rate  response 
to  a  step  aileron  input  is  shown  in 
Figs.  21(a)  and  21(b).  The  transfer- 
function  model  shown  in  the  dashed 
line  accurately  predicts  that  the 
response  will  be  predominantly  first 
order,  with  no  overshoot  or  oscilla¬ 
tory  tendency.  This  corroborates  the 
same  conclusion  made  earlier  on  Che 
basis  of  Che  small  and  1 4/6 1.1 

ratios.  The  sideslip  response  to  a 
step  rudder  input  is  shown  in 
Figs.  22(a)  and  22(b)  with  the  appro¬ 
priate  corrections  for  Che  nose-boom 
position.  The  transfer-function  model 
shown  in  the  dashed  curve  matches  the 
flight  data  fairly  well  for  the  first 
10  sec.  Thereafter,  the  responses 
diverge,  probably  due  Co  a  gust 
encounter  or  soite  other  aerodynamic 
interference  effect  at  Che  sideslip 
vane  during  Che  recovery  phase.  Once 
again,  the  character  of  Che  step 
response  is  well  predicted,  with  an 
accurate  modeling  of  the  rise  time, 
overshoot,  and  steady  state  sideslip. 

This  verification  study  shows 
the  utility  of  fairly  simple  transfer- 
function  models  in  predicting  the 
transient  responses  to  relatively 
large  and  varied  inputs.  High- 
frequency  and  mid-frequency  spectral 
response  matching  is  clearly  essential 
for  accurate  step  response  prediction, 
as  it  is  for  handling  qualities 
assessments.  Time  response  matching 
for  unstable  transfer  functions  is 
very  difficult  because  of  the  inabil¬ 
ity  Co  acquire  and  maintain  a  steady 
Initial  trim,  and  Che  tendency  toward 
rapidly  growing  transients  due  Co 
small  unmodeled  excitations. 


Figure  20.  Cosparison  of  aircraft  and  transfer- 
function  s»del  responaaB  to  a  step  elevator  input  in 
cruise,  (a)  Elevator  surface  deflection,  (b)  Fitch- 
rate  response,  (c)  Normal  (positive  downward)  accel¬ 
eration  response. 


-  FLIGHT  DATA 

- TRANSFER  FUNCTION  MODEL 


Figure  21.  Comparison  of  aircraft  and  transfer- 
function  ludel  responses  to  a  step  aileron  input  in 
cruise,  (a)  Aileron  aurface  deflection,  (b)  Noil- 
rate  response. 
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Figure  22.  Coaparlson  of  aircraft  and  transfer-function  model  responses  to  a  step  rudder  input 
in  cruise,  (a)  Rudder  surface  deflection,  (b)  Sideslip  response  (at  the  c.g.). 


S.  Conclusions 

Some  specific  conclusions  from  this  study  of  the  open-loop  dynamics  of  the  XV-15  tilt- 
rotor  aircraft  are; 

1)  The  response  characteristics  for  the  cruise  flight  condition  are  very  stable  and 
decoupled.  The  handling  quality  parameters  meet  or  exceed  Level  II >  Category  k,  requirements 
for  fixed-wing  military  aircraft. 

2)  Standard  lower-order  equivalent  Mdels  adequately  match  the  identified  cruise  dynam¬ 
ics  in  all  degrees-of-freedom  except  for  the  normal  acceleration  responses  to  elevator  deflec¬ 
tion.  The  misMtch  for  this  response  is  probably  due  to  the  unaK>deled  rotor  rpm  dynamics. 

3)  The  step  responses  of  the  Identified  transfer  functions  and  the  aircraft  generally 
match  very  well  for  both  the  hover  and  cruise  flight  conditions.  Time-domain  verification  of 
Che  unstable  attitude  dynamics  in  hover  is  very  difficult  because  small  errors  at  the  beginning 
of  the  run  cause  the  responses  of  the  model  and  the  aircraft  to  diverge  rapidly. 

The  frequency-dcauin  approach  has  proven  to  be  a  relatively  simple  and  accurate  means 
for  extracting  the  bare-airframe  dynamics  of  the  XV-IS.  The  utility  of  the  derived  lower  order 
transfer-function  andels  for  handling  qualities  and  control-system  studies  has  been  shown  in 
the  frequency  and  time  domains. 
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